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ABSTRACT 


A fomizentall circular cylinder, located near the floor of a two 
dimensional wave channel, was subjected to a train of gravity waves. The 
wave height and horizontal and vertical forces on the cylinder due to the 
incident waves were measured and presented in dimensionless form. Experi- 
mental values of the horizontal and vertical dimensionless force coeffi- 
cients are presented as functions of relative wave height and relative 
wave length. ‘The dimensionless force coefficients predicted by a modified 
Morrison's equation are compared with the experimental data. Experimental 
wave lengths of fran two to seventeen feet were investigated. 

The horizontal force coefficients were found to vary linearly with 
wave height over the range of wave lengths tested. The vertical force 
coefficients displayed regions in which the force was both inertia 


dominated and dominated by a lift force. 
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I. INTRODUCTION 


Exploration of the ocean floor has received increased interest in the 
past several decades, particularly the area of the continental shelf. A 
Significant amount of attention has been devoted to offshore oil explo- 
ration, to the recovery of petroleum and to deployment of undersea 
habitats. Often petroleum is located in remote areas of the world and 
oil production involves use of very large submerged oil storage tanks and 
associated pipelines. In other locations petroleum has been piped ashore, 
through the surf zone, to shore facilities. Recently, interest has been 
demonstrated [1] in the deployment of large submerged nuclear power 
plants. In the reference cited it is proposed that a nuclear plant be 
placed on the ocean floor in two hundred and fifty feet of water. 

As one might expect, this activity has generated considerable interest 
in the interaction of gravity waves with submerged objects. For example, 
the wave force exerted on vertical piling has been a topic of many investi- 
gations during the past twenty years [2], [3], [4]. In these studies the 
well known Morrison Equation [4] appears to be the most practical approach 
to the calculation of wave forces. This equation combines both a drag and 
inertia term, the values of which must be determined experimentally. At 
best, these two experimental coefficients are dependent upon the geometry 
of the pile and amplitude of the fluid motion, and it would appear that 
twenty years of research has only partially delineated this interdependence. 

In the case of large, deeply submerged objects the Morrison Equation 
is also applicable and generally yields valid results for the horizontal 
wave forces provided diffraction effects do not become appreciable. In 


this instance, the drag term may be discarded since the hydrodynamic force 





is primarily inertial in nature. However, it is still necessary to know 
the added mass coefficient as it depends on the body configuration and 
the effect of the rigid bottom surface. Moreover, if the object is 
bottom mounted the vertical fluid acceleration is zero, and accordingly, 
the Morrison Equation yields zero vertical force, a result that is 
obviously invalid. 

A diffraction theory, valid for the calculation of forces exerted on 
large submerged or semi-submerged bodies of arbitrary shape, which accounts 
for both the effect of the free surface and bottom, has been developed by 
Garrison and Chow [5]. Theoretical results were compared with corres- 
ponding experimental results for two different practical submerged tank 
configurations and good agreement was found. However, the method is valid 
for linear incident waves only; the nonlinear effects of large amplitude 
waves can not be determined. 

In the case of large objects, which are not so deeply submerged that 
the. free surface has a sizable effect, no theory is available which is 
valid for large values of wave height and few experimental results have 
been published. However, one study has been carried out by Johnson [6] 
for horizontal wave forces acting on a bottom-mounted, horizontal 
. Circular cylinder in long waves. Inasmuch as only long waves were con- 
Sidered, the results were dependent on the wave height and water depth 
only and no consideration was given to intermediate wave lengths. More- 
over, rather than correlate the forces with appropriate dimensionless 
parameters, Johnson used a regression analysis to'correlate the forces. 


Consequently, little light was shed on the basic mechanism involved. 





Schiller [7] measured wave forces acting on a submerged horizontal 
circular cylinder due to rather small amplitude waves. Over the range of 
wave heights considered in his tests, he found the amplitude of the hori- 
zontal force to vary linearly with the wave height. At small wave heights 
the vertical force was found to vary harmonically with the wave motion, 
the maximum amplitude of which increased linearly with the wave height. 

At larger values of wave height, the vertical force amplitude tended to 
increase in proportion to the wave height squared. This latter variation 
of the wave force is caused by the velocity squared term occuring in 
Bernoulli's Equation. 

The present investigation was an extension of Schiller's work to 
include large wave heights. However, the present investigation was limi- 
ted to the single configuration where the horizontal circular cylinder 
was placed on the bottom. A transition section was placed in the wave 
channel used by Schiller, at a location twenty feet from the wave 
generator for purposes of providing a natural transition from intermediate 
or deep water waves to shallow water waves of finite amplitude. The beach 
slope was adjusted for negligible reflection as discussed in the section 


on experimental apparatus. 
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II. THEORETICAL ANALYSIS 


The problem under consideration is depicted in Figure l. A train of 
regular waves of height H is considered to progress in the positive x 
direction in water of depth h. It is the present interest to determine 
the horizontal and vertical components of wave force acting on the 
horizontal circular cylinder in contact with the rigid bottom. 

The exact analytical solution to this problem is extremely difficult 
and, therefore, only an approximate analysis is attempted. However, it 
is first instructive to carry out a dimensional analysis of all of the 
pertinent parameters. It is known 4@ prtort that the maximum horizontal 
or vertical force per unit length of a cylinder is dependent upon the 
following variables: 


E 
x max y max 
== Clg 


R R 





= £(h,a,L,H,p 1S7H) (1) 


in which 


wave height 

cylinder length 

water depth 

cylinder radius 

wave length 

fluid density 
gravitational acceleration 
fluid viscosity. 


rave yee 
rbrk ut tow il 


It may be noted at this point that a relationship exists between the 
parameters associated with the incident wave (i.e., H, T, L and h, where 
T denotes the period). Consequently, only three of the incident wave 
parameters are needed in the dimensional analysis, i.e., either H, h and 


ion oH, h and T. Moreover, either of these three sets of parameters 


aie » 








FIGURE 1 DEFINITION SKETCH 
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campletely describes the incident wave regardless of whether or not it 
is of small amplitude. For small amplitude waves the well-known 
relationship . 


(2) 


exists between the parameters which describe the incident wave. For waves 
of finite amplitude the relationship is more complicated but nevertheless 
exists. 

A dimensional analysis of the physical paramaters indicated in 


Equation (1) gives the following dimensionless parameters. 





F F 
x Max y max _ ¢ aha a = u (3) 
pga22 pga72 L a 2a pyYgha2 


As pointed out by Garrison and Chow [5] the last temm on the right hand 
side of Equation (3) represents the ratio of the Froude number to 

Reynolds number and, accordingly, is an indicator of the ratio of the 
viscous to gravity forces. If this number is small, corresponding to 
large scale flows with small viscosity, it may be assumed that this 
parameter may not be important since the flow would be controlled primarily 


by gravitational forces. In this case Equation (3) may be written as 





F F 
oe eee Sta (4) 
pga22 pga22 i aoe za 
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However, as noted previously, the wave period may be used in place of the 


wave length, in which case Equation (4) would appear in alternate form as 





F F 2 
x max 4, ymax, jg h (5) 
pgaZe pga22 he ae 2a 


Equation (5) may also be written in a third form as 








F F 2 2 
ee a ee (6) 
pga22 pga22 h Hea 


In the experimental program an attempt is made at determining the corre- 
lation of the dimensionless parameters, on the right hand side of 
Equations (4), (5) and (6), with the force coefficients. 

The camplete analytical solution to the wave interaction problem is 
quite complex. However, by use of Morrison's Equation [4], an approxi- 
mate representation may be carried out. For this purpose, the velocity 
potential associated with the incident wave is written as 

p= Ho Cosme tye Castkx — ot) (7) 
296 «Cosh. ih 
where k = 27/L denotes the wave number and o = 21/T denotes the frequency. 
From Equation (7) the horizontal component of velocity and acceleration 


is obtained, respectively, by differentiation as, 


u = - # gk Cosh k(yth) chine = ele) (8) 
2o Cosh kh 

gat gy, COSh kth) cost - ot) (9) 
2 Cosh kh 
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The horizontal force acting on the cylinder is written, according to 
Morrison's Equation, as 


C 
ips a op 2ag u2 + (1 + om) Ta eu (10) 
2 


where Ca denotes the drag coefficient and om denotes the added mass 
coefficient. However, for cases where the particle motion is small in 
comparison to the cylinder diameter separation effects are small and it 
is possible to disregard the drag contribution to the total force and 


write Equation (10) as 


FY H Cos ( ot) 


= (fe + Cy Kk — —————— Gal} 
pga2Z 2 Cosh kh 





Using the definition of k and o, Equation (11) may be written in tems 


of the dimensionless parameters indicated in Equation (4) as 


Fy. max 


27a 7 
oga?2 (H/2a) 


L Cosh» ( 


= (1 +¢) 


n (12) 





27a ny 
io a 
The vertical component of the force could be expressed in a form 
similar to Equation (10), but the vertical camponent of acceleration on 
the bottom is zero. It is supposed, therefore, that the vertical 
component of force is associated with the lift force only. Accordingly 


the vertical force may be written as 


ene 
Fy sa C,, 2a2 (13) 


i 





in which C. denotes the lift coefficient and u the wave induced particle 
velocity on the bottom. Substituting for the velocity from Equation (8), 


Equation (13) may be written as 





II aon EE 
FY =2¢, 2ma\ {H_\ “ Sin*ot _ (14) 


pga22 L 2a Sinh 2kh 
Equations (12) and (14) contain two coefficients, namely the added mass 
coefficient, CA and the lift coefficient, C.. Values for those coeffi- 
Cients which account for both the bottom and free surface are not 
available in the literature. However, if the water is relatively deep 
the free surface effect should be small. Im this case the exact potential 


flow value for a circular cylinder in contact with the rigid bottom has 


been calculated by Garrison [8] as 


as = 2.29. 


For the same geometry, Daltman and Helfinstein [9] have calculated the 


lift coefficient by an approximate method. Their value is 


C= 4.48. 


These values of the two coefficients may be used in Equations (12) and 
(14) to obtain approximate expressions for the horizontal and vertical 
force coefficients. It should be noted, however, that the assumptions 
underlying these approximate relationships are met only when the scale 
of the waves, i.e., the wave length and water depth are large in compari- 
son to the cylinder diameter. Moreover, since linear wave theory has been 
used to determine the velocity and acceleration, the validity is also 


restricted to small amplitude waves. 
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III. DESCRIPTION OF APPARATUS AND EXPERIMENTAL PROCEDURE 


A. APPARATUS 

In order to measure the wave forces acting on the horizontal cylinder, 
experiments were carried out in a small wave channel. The basic wave 
channel, as shown in Figure 2, consisted of three sections; a deep water 
section, a shallow water section and a transition section, each fifteen 
inches wide with a rectangular cross section. The overall channel length 
- sixty-four feet. 

Three quarter inch exterior plywood sheets were used for the channel 
sides. Vertical rigidity was provided by bracing the plywood sheets with 
2 x 4's at various intervals, every four feet in the shallow water section, 
every foot in the deep water section and every two feet in transition 
section. Two 2 x 4's were attached to the top of the vertical studs and 
extended the length of the channel. 

Since the bottom of the deep water section was resting on a concrete 
foundation, only a single section of plywood was used for the channel 
floor. A double floor was used in the shallow and transition sections. 
This was done to provide extra strength. The double floor was constructed 
from two plywood sheets and two and a half inch separators. The sides of 
these sections were bolted together, through the spacing of the double 
bottom, with three-eighths threaded rod. Steel angle spacers were placed 
across the channel top to maintain the width dimension. 

Prior to assembly all wooden sections exposed to water were water- 
proofed. The channel interior was then given two coats of Morewear 
Vitri-Glaz 1320A. Dow Corning Sealant 780 was applied to all seams and 


joints. 
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1. Deep Water Section 


Several methods of increasing wave height and steepness are 
available. One of the most common methods is the generation of waves in 
a deep water regime and allow them to transition, by shoaling, toa 
shallow water regime. By the proper choice of depths and shoal slope, the 
transition can be such that there is no reflection from the shoal, and the 
resultant shallow water waves are in a near breaking condition. It was 
this method of increasing wave height that was chosen. As shown in 
Figure 2, the deep water section was Sixteen feet long and four feet deep. 
A wave generator was located two and three-quarters feet from the closed 
end of the channel. In order to generate regular waves in the channel a 
paddle type wave generator was used. An aluminum plate, adequately 
reinforced, was hinged at the wave channel floor and attached by a pin 
connection at the top to a driving rod, Figure 3. OQne-eighth inch teflon 
sheeting, four inches by fifty-one inches was located along the vertical 
sides of the plate to prevent leakage from one side to the other. In the 
space behind the wave paddle the wave motion was damped by use of baffle 
plates. A two horsepower variable speed drive, with an output speed 
range of twenty to one hundred and eighty revolutions per minute, was 
mounted on the top of the wave channel and fitted with a six and three- 
quarter inch radius face plate. The face plate was provided with ways 
so that the wave generator driving rod could be adjusted to any 
eccentricity between zero and eight inches, Figure 4. 

2. Transition Section 

As previously noted, the wave height was increased by shoaling 

deep water waves. To accomplish this a 14.14 foot long ramp was 


constructed and installed between the deep and shallow sections. This 


Ne) : 








FIGURE 3. WAVE GENERATOR PADDLE 
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provided for the transition from deep water to shallow water. In order 
to reduce the amount of reflection from the transition ramp, the slope 
was set at 1:6. The change in depth due to the shoal was two feet. 
3. Shallow Water Section 

The shallow water section had a depth of two feet and a length 
of forty-two feet. The test module was located twenty four feet from the 
transition; thus, allowing enough length for the waves to reach a fully 
developed state following the transition, before reaching the test section. 
At the end of the channel was located a variable slope beach for purposes 
of dissipating the wave motion. The beach consisted of metal shavings 
held between two sheets of perforated metal. These were separated by 
two inch wooden spacers. A solid sheet of one-eight inch aluminum was 
attached to the bottom of the beach by three and five-eighths inch 
separators. On the exposed surface of the beach was attached, parallel 
to the wave fronts, one by one inch aluminum angles. The maximum slope 
of the beach was 1:7. 

a. Test Module 

In order to provide a means of locating the test cylinder and 

measuring the wave forces acting on it, a separate test module was 
designed. The test module, consisting of the circular cylinder along with 
its own walls and floor, was constructed as an integral system in order 
that it may be removed from the channel. This concept allowed for the 
maintenance of close tolerances during manufacture and provided for easy 
access. ‘To accept the test module, the wave channel sides and floor were 


recessed and made of plexiglass. 
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The four inch diameter plexiglass cylinder was suspended 
between the walls of the test module by use of cantilever beams and 
adjusted to approximately one-sixteenth inch above the floor, Figure 5. 
It was known a priori that this sixteenth of an inch clearance would 
cause a high velocity jet under the cylinder due to wave action. To 
prevent this, a flexible plastic barrier was installed in the gap. The 
barrier was held in place by 'O' ring material pressed into 0.007 inch 
slot in the cylinder and the test module floor. 

In order to suspend the cylinder a small distance off the 
wave channel floor and at the same time measure forces exerted on it by 
surface waves, cantilever beam mounts were used as supports as shown in 
Figures 6 and 7. Bulkheads were fixed in the cylinder at approximately 
two inches from each end, and the fixed ends of the cantilever beams 
bolted on these bulkheads. The free ends of the two beams were allowed 
to protrude slightly beyond the end of the cylinder and into the plexi- 
glass test module walls. The free ends were supported by small self- 
aligning ball bearings pressed into the plexiglass. Both beams were 
fitted with strain gages and waterproofed using BLH Barrier 'C', 

Figure 8. One of the beams with its largest cross sectional dimension 
horizontal was used to measure the vertical force and the other with its 
largest cross sectional dimension vertical was used to measure the 
horizontal force. Calibration tests showed the cross coupling to be 
negligible. 

This design presented some unique problems. The cantilever 
beams had to be flexible enough to allow measurement of the forces and 


stiff enough so that the natural frequency was large in comparison to 
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HORIZONTAL FORCE CANTILEVER 


FIGURE 6. 


2D 








FIGURE 7. VERTICAL FORCE CANTILEVER 
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FIGURE 8. FORCE CANTILEVER 
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the excitation frequency. The characteristics of the strain gage bridge 
and amplifier determined the minimm strain necessary to produce a 
measurable output. The strain resulting on a cantilever from a given 


load is 


(15) 


where 


= strain 

= load 

= cantilever beam length 

modulus of elasticity 
cantilever beam base dimension 
cantilever beam height dimension 


Oe 6 
i 


The relationship for the spring constant of a cantilever beam is 


pts 
“ae 


K 





(16) 


By substituting Equation (15) into (16) a relationship for the spring 


constant, in terms of the strain and cantilever dimension is determined 


3 
Be i 
e°E*b*h, 


(17) 


It is apparent that the maximum stiffness and, consequently, the maximum 
natural frequency is obtained by making b, h and ec as small as possible. 


The minimum value of b was set by the width necessary to mount the strain 
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gages. The minimum valve of e was determined by the maximum sensitivity 
of the amplifier/recorder. The maximum stifMess was then achieved by 
making h as small as possible. However, Equation (15) shows that the 
beam length decreases with b and h as expressed by 

eEbh? 


4, ; 


6P 





Thus, the lower limit on h was finally determined by setting a reasonable 
minimum value for the beam length. 
In order to prevent the cylinder from rotating, an offset arm 
was connected from the horizontal force cantilever to the test module. 
A Farber bearing, AMS 1K7, was used as a wheel at the test module end of 
the arm to reduce friction. 
4, Wave Height Probe 

There are two types of probes in cammon usage for measuring wave 
height; these are commonly known as a capacitance and a resistance type 
probe. The resistance type utilizes ‘two vertical wires placed a small 
distance apart; the amount of wetting of the wires varies the resistance 
between them. This active element or "variable resistance" is wired 
into one leg of a Wheatstone bridge and the signal read out on an 
amplifier/recorder. Experience has shown, however, that the resistance 
type probe is extremely sensitive and well suited to measuring small 
amplitude waves but is rather nonlinear at larger amplitudes. 

The capacitance type wave height probe, although less sensitive, 
tends to be linear over large wave amplitudes and was, therefore, used 


as the primary method of sensing wave height in the present study. The 





probe consists essentially of a single insulated wire which acts as a 
variable capacitor; the wire acts as one plate and the water the second 
plate. As the water surface rises and falls, the capacitance varies 
accordingly. 

The probe used in the present study is shown in Figure 9. Located 
at the bottom of the foil-like support member was a three-eighths inch 
diameter acrylic rod. The sensing wire was terminated at the bottom in 
the acrylic rod and at the op at a nylon block. The tip of the acrylic 
rod was made removable in order to facilitate changing the sensing wire. 

A number of wire types were considered and tested prior to 
selecting number 30 A.W.G. wire with polythermaleze insulation. This 
Wire provided the capability of measuring wave heights to less than one | 
thirty-second of an inch with good accuracy. 

The electronic circuit utilized with this probe was obtained 
in schematic form from Hewlett Packard and is shown in Figure 10. A 


Hewlett Packard carrier amplifier 350 was used with the probe circuit. 


B. TEST PROCEDURE 

With the channel filled to the desired level, the test section and 
wave height probe were set in place. In each instance the probe was 
immersed to a depth of nine inches. The amplifiers were then adjusted 
to balance the bridge. 

With the use of a transverse mechanism the calibration of the wave 
probe was accomplished. The probe immersion was varied over a seven 
inch range in half-inch increments. The gage output was measured both 
on the Hewlett Packard recorder and one channel of a Brush recorder. 
This provided the initial calibration information. Prior to and after 


each set of runs the calibration was checked at several points. 
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FIGURE 10. WAVE HEIGHT CIRCUIT SCHEMATIC 
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The force cantilevers were calibrated by exerting a load on the 
cylinder by means of a series of weights and a wire pully arrangement as 
shown in Figure 11. The force was transmitted to the cylinder through a 
series of tapped holes located around the circumference at three, nine 
and twelve o'clock. Calibration information was obtained from the Bush 
recorder output by loading the cylinder in half-pound increments. As 
with the wave height probe, calibration was checked prior to and after 
each series of runs. 

With the calibration completed, the wave generator speed was adjusted 
for the desired wave length. The wave height probe was located an integer 
number of wave lengths fram the cylinder. Its output was connected to 
the horizontal grid of an oscilloscope. A second probe was then located 
directly over the cylinder and its output connected to the vertical 
oscilloscope grid. By observing the scope pattern the first probe could 
be adjusted to be in phase with the waves at the cylinder. The second 
probe was then removed and a data run commenced. During each run the 
wave generator drive eccentricity was varied from a half-inch to eight 


inches, or to a point where the waves broke in the channel. 
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FIGURE 11.CANTILEVER CALIBRATION ARRANGEMENT 
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IV. PRESENTATION OF RESULTS AND CONCLUSIONS 


A. EXPERIMENTAL RESULTS 

As discussed in Theoretical Analysis, non-dimensional wave forces can 
be represented as functions of relative water depth (h/a), relative wave 
length (2na/L) or period parameter (h/gT*) and relative height parameter 
(H/2a). In order to ascertain the effect of varying these parameters on 
the horizontal and vertical components of wave force a test program was 
established which, for a given run, held water depth and period constant 
while varying the wave height. Four series of runs were conducted at 
different water depths corresponding to the relative water depths (h/a) 
of 9.0, 7.0, 5.5 and 4.0. For each of these water depths, the relative 
wave length (2na/L) was varied over a range fran 0.06 to 0.48. 

Due to the difficulty in presenting all the wave force traces, a 
series of representative traces is shown in Figures 12, 13, 14 and 15. 
These traces represent three wave lengths from the three deeper water 
depths (h/a = 9.0, 7.0 and 5.5) and two wave lengths for the shallowest 
water depth (h/a = 4.0). For each of the eleven runs, four wave heights 
arranged in increasing order are shown. The uppermost trace represents 
wave height, the middle vertical force and the bottom trace is horizontal 
force. 

The wave height traces show that the waves were generally sinusoidal 
in nature. This was especially true at the shorter wave lengths. At the 
greater wave lengths the wave trace departed from the sinusoidal form as 
the wave height was increased. As can be seen in Figure 12, for the 
greater lengths, as the height was increased the wave was characterized by 


sharp peaks and long flat troughs. 
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Before discussing the wave force traces a brief discussion of the 
mechanism involved in the wave forces is in order. In Bernoulli's 
equation two contributions to the pressure acting on the surface of the 
cylinder are recognized, an unsteady contribution associated with the time 
derivative of the potential to the first power (vertical forces) and a 
tem associated with the fluid velocity squared. The velocity squared 
term tends to be symmetric about the y axis and, consequently, contributes 
little to the horizontal force. The inertia term is therefore dominant, 
in the horizontal forces, and increases linearly with wave height H/2a. 
Examination of the representative traces shows that the horizontal force 
varied nearly sinusoidally with wave height for all cases except when the 
wave length became very large. In the case of very long waves, the hori- 
zontal force peaks are separated by a null region which corresponded to 
the wave trough. The maximum value of the horizontal force coincided with 
the maximum water particle acceleration associated with the incident wave, 
Equation (9). 

For the vertical forces, there are two regimes, one in which the 
inertial forces dominate and a second in which the velocity squared term 
dominates. The inertial contribution to the pressure produces both posi- 
tive and negative values of vertical force. This term tends to be linear 
in wave height H/2a and have a maximum value of uplift when the cylinder 
is under the wave trough. The velocity squared term produces only upward 
forces due to the non-synmetrical flow about the cylinder. This component 
of vertical force is a function of the wave height squared (H/2a) 2 and 
peaks occur at the points of maximum wave induced velocity occurring at 
both crest and trough of the wave. Examination of the vertical force 


traces shows that for the longer waves the force was regular and hada 
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frequency twice that of the incident wave. The Terie frequency was 
caused by the velocity maximums at the crests and troughs. The maximm 
occurring under the wave crest in the longest waves was apparently 
attributable to the fact that the velocity under the crest of a shallow 
water wave is greater than at ae trough. Here the upward force was lift 
dominated. At intermediate wave lengths the reverse occurred. As the 
wave length was shortened the inertial effects gradually took effect and 
the maximum occurred at the trough where the two components were additive. 
Finally, at very short wave lengths the inertial force became dominant 
and the traces showed a sinusoidal behavior; the lift force associated with 
the velocity squared term in Bernoulli's equation became negligible. 

1. Horizontal Forces 

As noted previously, the horizontal force was nearly sinusoidal for 

the entire series of data runs. It is therefore possible to characterize 
the force by the maximum value. In order to examine the variation in the 
horizontal force with varying wave height, plots were constructed of the 
non-dimensional force coefficient fx (fx = I max! °98"*) versus the 
relative wave height H/2a. These wave force versus wave height plots are 
presented in Appendix B. It is apparent from the plots that horizontal 
force coefficient is characterized by a linear variation with wave height. 
Inasmuch as the variation of the force coefficient was linear with H/2a, 
for a given wave length, it is possible to characterize the horizontal 
force by the dimensionless parameter fx/(H/2a), the slope of the wave 
force-height plots. 

The parameter, ae ae, [pga* (H/2a)] is presented as a function of the 
relative wave length (2ra/L) in Figure 16. This figure also presents the 


result predicted by Morrison's analysis, Equation (12). For the larger 
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values of. relative wave length (21a/L) the experimental results and 
theory are in good agreement. This is generally as expected because for 
the larger values of 27a/L the waves are well represented by linear theory 
and the effect of the free surface is reduced due to the greater water 
depth. It may be recalled that these were the assumptions upon which the 
Morrison equation is based. However, for smaller values of 2na/L, the 
experimental results are consistently high. This is not unexpected as the 
waves have become shallow water waves and are not well represented by the 
linear wave theory. Also, in the case of the smaller depths the effect of 
the free surface becomes important. 

2. Vertical Forces 

As previously discussed, the vertical force may display two regimes, 
one inertia dominated and one dominated by the force associated with the 
velocity squared term in Bernoulli's equation. The assumption having been 
made in Equation (13) that the inertia forces would be negligible, the 
vertical forces would be a function of the velocity squared only and hence 
a function of (H/2a)2. It is difficult to characterize the force vari- 
ations with wave height by a single tem as two types of flow exist. 
However, having made the assumption that the force was a function of 
velocity squared, it was decided to use fy/ (H/2a)2 as the characteristic 
quantity. 

Values of Fy peep ga?2(H/2a) 2 were plotted against the relative wave 
length (27a/L) and presented in Figure 17. The values of this coefficient 
predicted by a lift force only are also shown in this figure. The agree- 
ment between the theory and experimental results appears to be good for 
the greater water depth (h/a = 9.0). However, as the water depth was 


decreased the agreement became poorer and the experimental data displayed 
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higher force values than predicted by Equation (14); apparently owing to 
the effect of the free surface and the non-linear effects of the incident 
wave. Also, as the wave length decreased, or the parameter 27a/L increased, 
it may be noted that the agreement is poor. This is apparently due to the 
fact that the inertia component becomes daminant and this contribution to 
the force was neglected in Equation (13). 

As waves are often characterized by period rather than wave length, 
Figure 18 presents the force coefficient as functions of a relative period 


ti/oT*).. 


B. CONCLUSIONS 

Fron the experimental results the following conclusions are considered 
warranted: 

1. The horizontal forces on the cylinder, resulting from the incident 
gravity waves, increased linearly with wave height. 

2. The vertical force showed two regimes, one inertial dominated and 
one velocity squared dominated. The inertial force regime occurred at 
shorter wave lengths, while the velocity squared term predcminated at the 
longer wave lengths. 

3. The maximum vertical force was always smaller than the maximm 
horizontal force and the downward vertical force was nearly zero. 

4, For greater water depths the theoretical predictions appeared to 
give valid results for both horizontal and vertical forces in the short 
wave length range. . 

5. Horizontal forces became independent of the parameter, 21a/L, for 


greater wave lengths. 
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TOTAL NUMBER OF RUNS IS 24. 
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TOTAL NUMBER OF RUNS IS 24. 
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TOTAL NUMBER OF RUNS IS 24. 
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NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 19. TOTAL 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 20. 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 21. 


WAVE PERIOD IS 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 22. 
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NUMBER OF RUNS IS 24. 


TOTAL 


0.678 SECONDS. 
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OF RUNS IS 24. 


NUMBER 


TOTAL 


0.655 SECONDS. 


DATA FOR RUN NUMBER 10- 
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DIMENSIONLESS PARAMETERS ARE 2PIA 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 25. 


WATER DEPTH IS 11.0 INCHES. 


WAVE LENGTH IS 13.46 FEET. 


2-549 SECONDS. 
DIMENSIGNLESS PARAMETERS ARE 2PIA 


WAVE PERIOD TS 


HO ett 
OMmost 
wor-0O 
Booe 
ee 8 @ 


Oooo 
ow te it 


DOODO0000 


@eeeeee# e@ 


OODO000°0 


mt OLAN + $0 
OOOKANNAAN 
eeeese eee 


QDOQO900000 


WW © LN OS 
AN STINOM oO 
eeeererecee 


COOoOOO00 


MORLINAMO Muy 
PAN AIO SF UNO 
eeeeete8ee®¢ 


OOOCO000°O 


ODOV00000 
evreeceecve 


ODOOCO000 


MODMAONAO 
OOAM + MOWO 
eeceecevere 


OCO3B0000°0 


MOYO OoMst 
MAM COOUTIO C 
eeeeveeee 


OOVOC HAHA 


aah ANAT D 
NOOSE OA 
®eee?eee? @ 


DOA AANA 


HNIC SLVOP- OO 
QMAIAIAIMAAYS 
LOD LA LALA EL LN 
MIAINAIONIOIAIN 


61 


XC 
‘N\ 








TCTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 27. 
1.436 SECONDS. 
DIMENSIONLESS PARAMETERS ARE 2PTA 
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WATER DEPTH IS 11-0 INCHES. 


WAVE LENGTH IS 7.09 FEET. 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 28. 


WAVE PERIOD IS 


WATER DEPTH IS 11-0 INCHES. 


WAVE LENGTH IS 5.36 FEET. 


1.150 SECONDS. 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 29. 


WATER DEPTH IS 11.0 INCHES. 


WAVE LENGTH IS 3.47 FEET. 


0.854 SECONDS. 


DIMENSIONLESS PARAMETERS ARE 2PIA 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 30. 


WATER DEPTH IS 11.0 INCHES. 


WAVE LENGTH IS 2.64 FEET. 


0.727 SECONDS. 


DIMENSIONLESS PARAMETERS ARE 2PIA 
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TOTAL NUMBER OF RUNS IS 24. 


RUN NUMBER 31. 


DATA FOR 


WATER DEPTH IS 11.0 INCHES. 


2.04 FEET. 


WAVE LENGTH IS 


0.634 SECONDS. 
DIMENSIONLESS PARAMETERS ARE 2PIA 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 35. 


WAVE LENGTH IS 7.30 FEET. WATER DEPTH IS 8.0 INCHES. 


1.658 SECONDS. 


DIMENSIONLESS PARAMETERS ARE 2PTA 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 36. 


Weve beNolhels Sa33 FEET. WATER DEPTH IS 8.0 INCHES. 


1.259 SECONDS. 
DIMENSTICNLESS PARAMETERS ARE 2PT 
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DATA FOR RUN NUMBER 37. TOTAL 
0.918 SECONDS. 
DIMENSTONLESS PARAMETERS ARE 2PIA 


WAVE PERTOD] Is 


NUMBER OF RUNS IS 24. 


3656 FEET. WATER DEPTH IS 8.0 INCHES. 


WAVE LENGTH IS 
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TOTAL NUMBER OF RUNS IS 24. 


DATA FOR RUN NUMBER 38. 


WAVE PERIOD IS 


WAVE LENGTH IS 2.50 FEET. WATER DEPTH IS 8.0 INCHES. 


0.724 SECONDS. 
DIMENSTONLESS PARAMETERS ARE 2PIA 
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